// Baccaro and Rei 2007 IO

// Load original dataset
clear
use "BR2007 IO Rep Data.dta"

// REPLICATION
summarize
display c(k)
mdesc
tabmiss
codebook country
codebook year

* TABLE 1
* column 1 pooled OLS
reg unr  lunr rir perdcpi tots ldprod ep ud_alt brr_alt  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_brr Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,
codebook year if e(sample)
codebook country if e(sample)

* column 2
xtgls unr  lunr rir perdcpi tots ldprod ep ud_alt brr_alt  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_brr Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,p(h)corr(ar1)
codebook year if e(sample)
codebook country if e(sample)

* column 3
xtgls unr  lunr rir perdcpi tots ldprod ep ud_alt brr_alt  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_brr Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,p(h)corr(psar1)
codebook year if e(sample)
codebook country if e(sample)

*column 4
xtgls unr  lunr rir perdcpi tots ldprod ep ud_alt benOECD tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep Dev_k_benoecd Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,p(h)corr(psar1)
codebook year if e(sample)
codebook country if e(sample)

* column 5
xtpcse unr  lunr rir perdcpi tots ldprod ep ud_alt brr_alt  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_brr Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,
codebook year if e(sample)
codebook country if e(sample)

*column 6
xtpcse unr  lunr rir perdcpi tots ldprod ep ud_alt benOECD  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_benoecd Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,
codebook year if e(sample)
codebook country if e(sample)

// Could not replication columns 7 and 8 due to absence of data

// IMPUTATION
// m = 5
clear
use "BR2007 IO Imp Data.dta"
summarize

mi import flong, m(imp) id(CNTRY year) imp(epop-Dev_ken_cbi benOECD-Dev_k_benoecd)

gen cntry = CNTRY

* TABLE 1
* column 1 pooled OLS
mi estimate, post: reg unr  lunr rir perdcpi tots ldprod ep ud_alt brr_alt  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_brr Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,
* column 2
mi estimate, cmdok post: xtgls unr  lunr rir perdcpi tots ldprod ep ud_alt brr_alt  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_brr Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,p(h)corr(ar1)
* column 3
mi estimate, cmdok post: xtgls unr  lunr rir perdcpi tots ldprod ep ud_alt brr_alt  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_brr Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,p(h)corr(psar1)
*column 4
mi estimate, cmdok post: xtgls unr  lunr rir perdcpi tots ldprod ep ud_alt benOECD tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep Dev_k_benoecd Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,p(h)corr(psar1)
* column 5
mi estimate, cmdok post: xtpcse unr  lunr rir perdcpi tots ldprod ep ud_alt brr_alt  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_brr Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,
*column 6
mi estimate, cmdok post: xtpcse unr  lunr rir perdcpi tots ldprod ep ud_alt benOECD  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_benoecd Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,

// m = 30
clear
use "BR2007 IO Imp Data 30.dta"
summarize

mi import flong, m(imp) id(CNTRY year) imp(epop-Dev_ken_cbi benOECD-Dev_k_benoecd)

gen cntry = CNTRY

* TABLE 1
* column 1 pooled OLS
mi estimate, post: reg unr  lunr rir perdcpi tots ldprod ep ud_alt brr_alt  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_brr Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,
* column 2
mi estimate, cmdok post: xtgls unr  lunr rir perdcpi tots ldprod ep ud_alt brr_alt  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_brr Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,p(h)corr(ar1)
* column 3
mi estimate, cmdok post: xtgls unr  lunr rir perdcpi tots ldprod ep ud_alt brr_alt  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_brr Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,p(h)corr(psar1)
*column 4
mi estimate, cmdok post: xtgls unr  lunr rir perdcpi tots ldprod ep ud_alt benOECD tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep Dev_k_benoecd Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,p(h)corr(psar1)
* column 5
mi estimate, cmdok post: xtpcse unr  lunr rir perdcpi tots ldprod ep ud_alt brr_alt  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_brr Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,
*column 6
mi estimate, cmdok post: xtpcse unr  lunr rir perdcpi tots ldprod ep ud_alt benOECD  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_benoecd Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,

// m = hd
clear
use "BR2007 IO Imp Data HD.dta"
summarize

mi import flong, m(imp) id(CNTRY year) imp(epop-Dev_ken_cbi benOECD-Dev_k_benoecd)

gen cntry = CNTRY

* TABLE 1
* column 1 pooled OLS
mi estimate, post: reg unr  lunr rir perdcpi tots ldprod ep ud_alt brr_alt  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_brr Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,
* column 2
mi estimate, cmdok post: xtgls unr  lunr rir perdcpi tots ldprod ep ud_alt brr_alt  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_brr Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,p(h)corr(ar1)
* column 3
mi estimate, cmdok post: xtgls unr  lunr rir perdcpi tots ldprod ep ud_alt brr_alt  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_brr Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,p(h)corr(psar1)
*column 4
mi estimate, cmdok post: xtgls unr  lunr rir perdcpi tots ldprod ep ud_alt benOECD tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep Dev_k_benoecd Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,p(h)corr(psar1)
* column 5
mi estimate, cmdok post: xtpcse unr  lunr rir perdcpi tots ldprod ep ud_alt brr_alt  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_brr Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,
*column 6
mi estimate, cmdok post: xtpcse unr  lunr rir perdcpi tots ldprod ep ud_alt benOECD  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_benoecd Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,

// m = mice
clear
use "BR2007 IO Imp Data MICE.dta"
summarize

mi import flong, m(imp) id(CNTRY year) imp(epop-Dev_ken_cbi benOECD-Dev_k_benoecd)

gen cntry = CNTRY

* TABLE 1
* column 1 pooled OLS
mi estimate, post: reg unr  lunr rir perdcpi tots ldprod ep ud_alt brr_alt  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_brr Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,
* column 2
mi xtset cntry year
mi estimate, cmdok post: xtgls unr  lunr rir perdcpi tots ldprod ep ud_alt brr_alt  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_brr Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,p(h)corr(ar1)
* column 3
mi estimate, cmdok post: xtgls unr  lunr rir perdcpi tots ldprod ep ud_alt brr_alt  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_brr Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,p(h)corr(psar1)
*column 4
mi estimate, cmdok post: xtgls unr  lunr rir perdcpi tots ldprod ep ud_alt benOECD tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep Dev_k_benoecd Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,p(h)corr(psar1)
* column 5
mi estimate, cmdok post: xtpcse unr  lunr rir perdcpi tots ldprod ep ud_alt brr_alt  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_brr Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,
*column 6
mi estimate, cmdok post: xtpcse unr  lunr rir perdcpi tots ldprod ep ud_alt benOECD  tw_alt2 cbi kenwcoor Dev_K_UD Dev_ken_tw  Dev_k_ep  Dev_k_benoecd Dev_ken_cbi ce2-ce20 ye2-ye39  if CNTRY!=15 & CNTRY!=16,

